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Abstract 



We recently derived, using diagrammatic methods, the leading-order hard 
photon emission rate in ultra-relativistic plasmas. This requires a correct 
treatment of multiple scattering effects which limit the coherence length of 
emitted radiation (the Landau-Pomeranchuk-Migdal effect). In this paper, we 
provide a more physical derivation of this result, and extend the treatment to 
the case of gluon radiation. 



I. INTRODUCTION 



The rate of photon emission in a high-temperature QCD plasma is a problem of some 
theoretical interest |l due in part to the hope that hard photon emission will be a 
useful diagnostic probe of heavy ion collisions |I5|,I6|. Understanding the analogous rate of 



gluon emission in a hot QCD plasma is required for computing thermalization and transport 
processes 



17-24 



We will make the simplifying assumption that the temperature T is so large that the 
running strong coupling constant g s (T) can be treated as small, and focus on the leading- 
order (in a s ) behavior of the photon or gluon emission rate.[] Evaluating just the leading 
weak coupling behavior, even in this asymptotically high temperature regime, is non-trivial. 
For the photon emission rate, one might expect that it would be sufficient to compute just 
the diagrams shown in Fig. [I], which describe lowest order two-to-two particle scattering 
processes. But it turns out that naive perturbation theory (augmented by the standard 
resummation of hard thermal loops) does not suffice to calculate the leading-order hard 
photon emission rate. Fig. ^| shows bremsstrahlung and pair production processes which 
also contribute to the on-shell photon emission rate at leading order. The essential problem 



■"■In other words, we will neglect all sub-leading corrections to the emission rate suppressed by additional 
powers of a s , as well as power corrections suppressed by Aqqu/T or m q /T. 
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FIG. 1. Two-to-two particle processes contributing to the leading order photon emission rate. Time 
may be viewed as running from left to right. 




FIG. 2. Brcmsstrahlung and pair production contributions to photon emission. The bottom line in 
each diagram can represent either a quark or a gluon. 




FIG. 3. An interference term, involving amplitudes for photon emission before and after multiple 
scattering events, which contributes to the leading order emission rate. 



is that the internal time scale associated with these processes is comparable to the mean 
free time for soft scattering with other particles in the plasma. Consequently, even at 
leading order, the photon emission rate is sensitive to processes involving multiple scatterings 
occurring during the emission process. This is known as the Landau-Pomeranchuk-Migdal 
(LPM) effect. Diagrammatically, it manifests as the presence of interference terms involving 
multiple collisions, such as depicted in Fig. [3], which parametrically are equally as important 
as those of Fig. [|. In the nearly-collinear limit, the extra explicit factors of g s turn out to 
be canceled by a combination of parametrically large enhancements from the internal quark 
propagators and soft exchanged gluons. 

In Ref. |p?5|| , we showed how to account for the LPM effect in photon emission by identi- 
fying and summing the appropriate infinite class of diagrams which contribute to the leading 
order emission rate, and in Ref. [[2(J we solved the resulting integral equations numerically. 
Our analysis in Ref. |25j involved quite detailed power counting of diagrams and relied upon 
some rather unintuitive relationships found by Wang and Heinz [27] between real-time ther- 
mal 4-point Green's functions in the r/ a (Keldysh) formalism. One goal of the present article 
is to reproduce our previous results in a manner that more clearly highlights the essential 
physics of the result. We will embrace physical arguments wherever possible, leaving to our 
previous work [^5| the more technical justification of our results. Our second goal will be to 
extend the treatment to the case of gluon emission. 

We are hardly the first people to discuss applications of the LPM effect, either qual- 
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itatively or quantitatively |p!3| , p!8H2^j2^ -p3| .p| However, discussions prior to our work f25 



have almost always focused on systems where the scatterers with which the emitting parti- 
cle interacts are static. (The scatterers are represented by the bottom lines in Figs. |2| and 
H) But typical scatterings in an ultra-relativistic plasma involve interactions with excita- 
tions which are themselves moving at nearly the speed of light. They produce dynamically 
screened color electric and magnetic fields which form a fluctuating background field in 
which bremsstrahlung, pair annihilation, and the LPM effect take place. 

Our discussion will be largely self-contained and will not directly rely on previous treat- 
ments of the LPM effect in other contexts. One could, presumably, explicitly mimic previous 
discussions such as the seminal 1955 analysis of Migdal p9| , |30| , and suitably generalize the 
treatment to non-static scattering. However, it is rather challenging to follow in detail all 
of the assumptions and approximations made in Migdal's quantum mechanical treatment. 
We believe that readers interested in understanding the effect can benefit from several dif- 
ferent formulations. Moreover, our discussion of scales and approximations will be tailored 
to the particular problem at hand, namely hard photon (and later gluon) emission in ultra- 
relativistic plasmas. 

We will focus on the contributions of bremsstrahlung and pair annihilation to the rate 
of photon (or gluon) emission. We will not explicitly discuss contributions to the emission 
rate from the 2^2 processes of Fig. [I], which are of the same order but do not require a 
treatment of the LPM effect. These 2^2 contributions to the photon emission rate are 
calculated in Ref. P, p!0|p6| f| Diagrammatically, the bremsstrahlung and pair annihilation 
processes of Fig. |^ look like they are suppressed by an explicit factor of gl compared to the 
2 «-> 2 processes of Fig. [I]. However, the diagrams of Fig. |2| have a collinear enhancement, 
associated with emission of the photon. This enhancement is produced by a would-be 
collinear singularity which is cut off by the effective thermal mass (of order g s T) and width 
of the fermions, yielding a parametrically large near on-shell enhancement which causes 
these processes to be the same order (up to a logarithm of g s ) as those of Fig. [I]. 

We will use the following nomenclature in our discussion. The particles corresponding 
to the lower lines in Figs. || or § will be called the scatterers. The particle which emits 
the photon in bremsstrahlung, or the particle/anti-particle pair which annihilate in pair 
annihilation, will be called the (photon) emitters. Each separate gluon exchange in Figs. |2| 
and [3] will be referred to as a "scattering" of the emitters.^ 

Before we begin, we first summarize the basic scales associated with bremsstrahlung and 
pair production of hard photons via the processes of Fig. |2|. For simplicity, we restrict atten- 
tion to photons with momenta parametrically of order T, which is the typical momentum of 



2 For a review of many aspects of the LPM effect, with an emphasis on the theory and experiment of 
high-energy collisions of electrons with atomic matter, see Ref. [ p4| . 

3 Only the k > T limit is addressed in Refs. 

4 When distinguishing emitters from scatterers, one need not worry about the fundamental indistinguisha- 
bility of quarks in, say, quark-quark scattering, because the scattering is dominated by small angle collisions, 
and the photon is emitted nearly collinear with the emitter(s). So, at leading order, emitters and scatterers 
are distinguishable simply by whether or not their direction of motion is nearly aligned with the emitted 
photon. 
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FIG. 4. Orders of magnitude of various momentum, distance, and angular scales associated with 
bremsstrahlung of a photon with momentum of order T. g stands for the strong coupling g s . 



particles in the ultra-relativistic plasma. We will throughout treat electromagnetic couplings 
as small compared to the QCD coupling, and we will ignore thermal effects on the propa- 
gation of the photon (which are proportional to a EM ). The basic kinematics, summarized 
pictorially in Fig. || for bremsstrahlung, turn out to be as follows. 

• The typical momenta of the emitters and scatterers is order T. 

• The typical momentum transfer q of an exchanged gluon responsible for scattering 
during the emission of the photon is of order g s T, which is also the order of the inverse 
Debye screening length for color. The angle of deflection in such a collision is of order 
this momentum transfer divided by the emitter momentum of order T, hence 9 ~ g s . 

• Processes contributing to the leading-order emission rate are dominated by nearly 
collinear emission of the photon. The corresponding internal quark lines in Fig. |2| 
are nearly on-shell, with energies off-shell by an amount 5E of order g^T. Fourier 
transforming, this implies that these processes have time durations of order (g^T)^ 1 . 
This is known as the formation time of the photon. 

• The formation time scale (g^T)' 1 is also the order of the mean free time between 
collisions with momentum transfers of order g s T. This is why multiple collisions cannot 
be treated independently and interferences such as Fig. || must be included at leading 
order. 

• The typical angle between the directions of the photon and the emitter(s) (be they 
initial state or final state emitters) is 9 ~ g s . 

• The typical angle between the directions of the emitter and a scatterer is 0(1). 



A brief review of how to obtain these scales is given in Sec. LB of Ref. p5 |. Here, we will 
take them as our starting point and proceed to discuss how to sum up the effects of multiple 
collisions. 

Certain physical processes in high temperature gauge theories (such as the rate of baryon 
number violation in hot electroweak theory) are sensitive to "ultra-soft" collisions which are 
mediated by low-frequency non-Abelian magnetic fluctuations with momentum transfers of 
order g%T. The dynamics of such ultra-soft collisions are intrinsically non-perturbative. One 
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may check a posteriori, by inspection of the final answers we will produce, that the leading- 
order dynamics relevant for hard photon (or gluon) production is not sensitive to ultra- 
soft collisions. Alternatively, the reader may find detailed qualitative and diagrammatic 
discussions of this point in Ref. |25 |. 

For simplicity, we will restrict attention to the case of zero chemical potential. The 
generalization of results for photon emission to non-zero chemical potential may be found in 
Ref. |E5|. In the next section, we will mainly focus on bremsstrahlung and reformulate the 
problem as that of bremsstrahlung from an emitter that is propagating through a classical 
random background color gauge field. We will show how simple physical considerations of 
localization of particles in both space and momentum lead to a description of the LPM effect 
in bremsstrahlung as an infinite sum of ladder diagrams. In section |ITI|, we will discuss how 
to convert the required sum of ladder diagrams into a Boltzmann-like integral equation — a 
task previously carried out in Ref. J25|, but which we will repeat using the formalism of this 
paper, in a way that more clearly displays which portions of the relevant physics depend on 
the nature of the emitter (e.g. an actual Dirac quark, a fictitious scalar quark, or something 
else) and which physics does not. In section |Vj, we then analyze the closely related process 
of pair production, and we combine these results for photon emission in section |Vj Finally, 
in section [VI] we discuss the generalization to gluon emission. 



II. RANDOM BACKGROUNDS AND ORDERING OF INTERACTIONS 

The differential photon emission rate (per unit volume), at leading order in a EM , is given 
by the well-known relation 

dT " = (2^k\ £ £ W (k) * ^ (k) W ^ {K) ' (2 ' 1} 

where K = (k°, k) = (|k|, k) denotes the null photon 4-momentum, the e's represent a basis 
of transverse polarizations for the photon, and W^lK) is the Wightman electromagnetic 
current- current correlator, 

W, U (K) = J d'x e~ lKx (j,(0)Ux)) . (2.2) 

We use a metric with signature ( — h++). Because the correlator is contracted with photon 
polarization vectors in ( [2.1|) , one need only consider the case where fi and v are spatial 
indices in what follows. As always, (• ■ ■) denotes an expectation value in whatever density 
matrix V is of interest, which in our case is a thermal ensemble describing the equilibrium 
plasma. If one inserts a complete set of multi-particle states |/) between the currents, and 
works in a basis |i) where V is diagonal, then the correlator can also be written as 

W, V {K) = J d'xe-^Y.^ (*L7m(0)I/> C/U(*)K>- ( 2 - 3 ) 

i 

It will be useful to remember that this corresponds to scattering from an initial state i to a 
final state / plus the emitted photon. Moreover, the integrand in (|2.3|) can be interpreted 
as representing interference between photon emission at the space-time points and x. 
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FIG. 5. The perturbative bremsstrahlung and annihilation processes of Fig. with the soft gluon fields 
now interpreted as classical background fields. 



+ 




FIG. 6. Bremsstrahlung and annihilation processes including the multiple interactions with a back- 
ground field, which lead to the LPM effect. 



Our first approximation will be to treat the soft (q ~ gT) gluons of the problem as 
a random, classical non-Abelian background field A^x), through which the particles that 
bremsstrahlung or annihilate propagate. The lowest-order processes of Fig. ^ are then re- 
placed by Fig. H and our task will be to sum up multiple gluon interactions such as depicted 
in Fig. ^[ After computing rates, we will appropriately average over this random classical 
gauge field. Such averaging will be denoted by ((• • •)). By translation invariance, the variance 
of the background field must have the formQ 

({A»{Q)*A»{Q'))) = P7(Q) {2^\Q-Q') ) (2.4) 

where p^ (Q) is the spectral density of thermal gauge field fluctuations. We may treat the 
statistical distribution of the background gauge field as Gaussian, dropping higher-order 
connected correlations. This is equivalent to neglecting the nonlinear self-interaction of the 
background and neglecting the interactions between different scatterers. This approximation 
is valid, at leading order, everywhere except in the deep infrared (Q ~ g 2 T). Neglecting 
non-Gaussian correlations in the background gauge field would not be allowable if the pro- 
cesses under consideration were sensitive (at leading order) to ultra-soft g 2 T fluctuations. 
Fortunately, this is not the case 



25p. 



Physically, these soft gauge fields are created by other charge carriers [with typical 0(T) 
momenta] passing randomly through the plasma. The statistics of the soft gauge fields 
generated by these particles can be described using (i) the fluctuation-dissipation theorem, 
and (ii) the standard hard thermal loop (HTL) approximation for the retarded self-energy 
n(Q) of soft (Q <C T) gauge fields. Namely, the spectral density Pa(Q) = \ \p I a(Q)\ \ ma Y be 
taken to be 



p A (Q)^2[ nb (g u ) + l]ImC L (Q) 



2[n b (g°) + 1] Im 



(2.5) 



5 A translationally invariant choice of gauge fixing, such as Lorentz or Coulomb gauge, is tacitly assumed. 
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where n\,((J) = l/(e^ a '— 1) is the Bose distribution function, G^ t (Q) is the free retarded 
gauge field propagator, and G^^(Q) is the retarded propagator with hard thermal loops 
resummed, which account for Debye screening and Landau damping of the color fields. We 
will review the specific formulas later, in section [V]. Because the background fields of interest 
are soft, one may replace the [ri B (<?°)+l] statistical factor by its low frequency limit, namely 
T/q°. 

The presence of correlations in the background field with non-trivial frequency (and 
wave-number) dependence is the most significant difference between our problem of hard 
photon emission in a relativistic plasma, and the original setting for the LPM effect of hard 
photon emission in a static medium of random Coulomb scatterers (nuclei). In order to be 
able to make contact later between our results and those of Migdal, we will treat the spectral 
density Pa{Q) as arbitrary, and not use the specific form ( |2.5| ) until the very end. 



A. Particle propagation in a random field 



Our next task is to rewrite the current correlator ( |2.2| ) in terms of propagation amplitudes 
of single particle states in the random background. We will focus on the bremsstrahlung 
process and come back to pair annihilation in section [TV]. Bremsstrahlung corresponds to a 
contribution to the current correlator ( [2.2|) of the formF] 

W]ZT{K) = J d'xe-^ J^n( Pi )[l±n( Pf )] i( Pi \j,(0)\Pf) (p f \jv(x)\ Pi )\ , (2.6) 

where |pj) and |p/) represent one-particle states of the emitting charge carrier, and n is the 
corresponding equilibrium distribution function [either n^(p) or nt(p)]. The contributions 
( |2.6| ) to the current correlator have to be summed over all types (including anti-particles) and 
spins of hard charged quasi-particles. The momenta Pi and pf can be regarded as the initial 
and final momenta of the charge carrier, and n(pi) and l±n(pf) are the initial and final state 
statistical factors. There is an approximation hiding here because momentum eigenstates 
|p) are not exact energy eigenstates in the presence of the soft background field. The density 
matrix is not therefore precisely n(p). Over the 0(l/g 2 T) formation time relevant to our 
problem, soft momenta transfers change p by order gT <C p. But the resulting error in the 
use of n(p) is then sub- leading in g and so can be ignored in a leading-order analysis. 

Making the time evolution explicit, and inserting complete sets of intermediate states, 
the statistically averaged matrix elements appearing in (|2.6|) can be rewritten as 



r pV ( (Vi\m\Pf) (Vf\U(0,x°)\p' f ) (p' / |j,(x)| P :)(p' t |f/(a; ,0)|p,^ 



pV <P/b'„(0)|p*r <P/|j*(x)|p£> ((p' f \U(x°,0)\p f y < P ;|f/Or ,0)|p 4 >) , (2.7) 



3 For any three-momentum p, J means J e? 3 p/ (2tt) 3 , while for a four-vector Q, J Q will denote J d 4 Q/ (27r) 4 . 
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FIG. 7. Two diagrams, time ordered from left to right, whose interference contributes to the rate of 
bremsstrahlung. The first diagram represents photon emission at time zero, and the second at time x°, and 
in both cases the diagrams show the evolution between these two times. 




FIG. 8. A time-ordered Z diagram. 




FIG. 9. A single diagram depicting the interference of the two diagrams of Fig. M. The interactions 
along top and bottom lines arc independently time ordered from left to right. 

where U(t',t) is the time evolution operator in the background field. As usual, the spatial 
Fourier transform in ( ^.6| ) combines with translation invariance (of the statistically averaged 
matrix elements) to enforce total momentum conservation: 

k = p. J -p / = P ;-p / / . (2.8) 

Diagrammatically, this contribution represents the interference in the evolution of the 
particle from time to x° depending on whether the photon is emitted at time or at time 
x°, as depicted in Figs. 0a and b. These diagrams can be considered as time ordered (with 
time running from left to right) because each individual Q ~ gT momentum transfer is 
not enough to create or destroy a (nearly on-shell) particle/anti-particle pair. That is, in 
a time-ordered Z contribution like Fig. |8|, the three-particle intermediate state would have 
to be so far off shell that its contribution is suppressed compared to the time ordering of 
Fig. [7|a. It is convenient to put the interference of the evolutions of Figs. [7| together into 
the single diagram depicted by Fig. The top line (or "rail") represents (p^|Z7 0)|pj) 
and the bottom line (rail) the complex conjugate of (pj|£/(x°, 0) |p/)- This looks just like a 
Feynman diagram for the current correlator, except with the added interpretation that each 
rail of the diagram can be considered as time ordered from left to right. 



S 
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FIG. 10. The emitter's trajectory, thought of as approximately classical, and the ordering in time and 
space this gives to its local interactions with scatterers (dashed lines). 



B. Ordering of interactions 

The next step is to understand the effect of averaging over the background gauge field. 
A quick way to visualize the dominant correlations is to think of the photon and emitter 
qualitatively as approximately classical particles, having both approximately well-defined 
position and momentum, obeying roughly 

x(t) ~ x(0) + v p t. (2.9) 



Indeed, the original, qualitative derivation of the LPM effect was purely classical p8|.f| 



(See section 1 of Ref. |24| for a concise review.) More formally, instead of using a basis 
of momentum or position states, one may use an (over-complete) basis of Gaussian wave 
packets with width Ax ~ (Ap) _1 ~ (gT)^ 1 . For momenta of order T, the widths will remain 
of this order over the time scale l/g 2 T relevant to our problem.^] The spatial uncertainty 
Ax in the position of the particles is then always small compared to l/g 2 T, which is the 
mean free path for the soft gluon interactions with Q ~ gT which dominate our problem. 
This implies that interactions with scatterers must occur in a definite order, regardless of 
when the bremsstrahlung photon is emitted. 

Fig. [10] illustrates the basic point. It shows the emitter propagating along the +z axis 
as a (nearly) localized classical particle. The dashed lines depict the small subset of other 
hard particles in the plasma which happen to pass close to the emitter during its flight 
and happen to interact with it when they do. Because of the localization of the emitter in 
space and time, and because of the local nature of the soft interactions (set by the Debye 
screening length 1/gT), the interactions (for the particular background of hard particles 



7 Classical results become precise quantitatively in the limit that the photon energy is small compared to 
the emitter energy. 

8 In more detail, the transverse momenta in the wave packet will be order Ap± ~ Ap ~ gT, both from 
the initial width as well as momentum transfer from scattering. The transverse spatial spread of the wave 
packet over time r ~ 1/ 'g 2 T will then be rAp±/p ~ 1/ 'gT ', since p ~ T. There is also a longitudinal spread 
in the wave packet due to the effective thermal mass ~ gT of hard quarks, which gives a dispersion in 
velocities of order Aw ~ A(p/ y/p 2 + m^J ~ m^Ap/p 3 ~ g 3 . The longitudinal spread of the wave packet 
due to this dispersion is then of order r Av ~ g/T. 
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FIG. 11. Two amplitudes which do not interfere (at leading order) because of mismatched order of 
encountering scatterers. 
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(a) fb) fc) 

FIG. 12. Amplitudes for which the time and space ordering of interactions is consistent. 




FIG. 13. The product of amplitudes, as in Fig. H, now averaged over random background fields. The 
dashed lines denote the independent correlations (\AA)) of the background gauge fields. The double lines 
represent the resummed propagators of Fig. H4. 



shown) must happen in the order numbered in the diagram, regardless of when the collinear 
photon is emitted. That is, there can be no interference (at leading order in g) between 
the diagrams of Fig. |TT] if the numbers 1 and 2 denote the fields generated by two different 
scatterers. However, the diagrams of Fig. fl2| a and b could interfere. So could the diagrams 
of Fig. |T2"| a and c, except that this particular interference vanishes when averaged over the 
randomness of the scatterers: the unmatched scatterer 4 could be a particle or anti-particle 
(for example), and a single insertion of its field in the product of amplitudes averages to 
zero by charge conjugation invariance. So, after averaging, non-zero contributions only arise 
from the cases where interactions with a given scatterer happen to appear twice in Fig. |9|, 
with the interactions of the top and bottom rails consistently ordered in time and space. We 
can rewrite the random average of these figures in the form of Figs. |13] and [14|, where the 
dashed lines indicate that the corresponding pair of A's have been replaced by the correlation 
i^A^Ay)). In our problem, this correlator is given by (|2.4j) , but the argument we have used to 
arrive at Fig. |1^ would also apply to static random scatterers, such as considered by Migdal, 
whenever one has a similar hierarchy of scales. 

We should mention in passing that the condition that the entire wave packet remain well- 
localized, in all three dimensions, on a scale parametrically small compared to the mean 
separation between scattering events, is much stronger than necessary to argue that the 
interaction must be ordered in time (at leading order). Interactions are effectively ordered 
in time, and ladder diagrams dominate, in numerous other applications including the classic 
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FIG. 14. The resummation of self-energy insertions into the propagator. 



case of scattering from point-like static random impurities.^ 

The self-energy resummation on the top and bottom rails, involving insertions of the 



background field correlator as shown in Fig. [TJ], represents the inclusion of the thermal 
width T in the propagator of the emitting particle. In hot gauge theories, the thermal 
width of quasiparticles is dominated by soft scattering and is actually an infrared divergent 
quantity in perturbation theory. However, these divergences (which represent sensitivity to 
ultra-soft fluctuations) will turn out to cancel in our final result when combined with the 
correlations that connect the top and bottom rails. Our final formulas will be well behaved 
in the infrared. 

Before moving on to the evaluation of the ladder diagrams of Fig. [13], we may use the 
fuzzy classical particle picture to understand why the LPM effect depends crucially on the 
fact that the photon and the emitting particle move at or very close to the speed of light. If 
a particle of momentum p moving in the +z direction emits a collinear photon of momentum 
k at time and position zero, then the subsequent positions of the particle and the photon at 
time t' will be 

•^photon — Vjt , ^particle — Vp—k^ i (2-10) 

where t> 7 is the velocity of the photon and v p the velocity of the (quasi) particle with 
momentum p. If, on the other hand, the same particle doesn't emit the photon until time t, 
the subsequent positions will be 

^photon ~V p t + V^ (t' - t) , ^particle ^ V p t + V p _ k (t' - t) . (2.11) 

These two processes can interfere significantly only if the photon and particle trajectories 
are the same, up to the fuzziness Ax. That requires 

Vryt ~ Vpt ~ Vp-kt, (2-12) 

which is natural (that is, not suppressed by additional phase space restrictions) only if the 
particles are all ultrarelativistic. 



9 In this case, a localized wave packet impinging on a single scattering center will produce a scattered wave 
which resembles an outgoing spherical shell, with no localization in direction whatsoever. But as long as the 
temporal duration of the scattered wave (as it moves past a fixed location) is small compared to the mean 
time between collisions, multiply-scattered waves produced by scattering, in different orders, off a given set 
of scattering centers, will have negligible overlap and hence not interfere. In particle-oriented language, all 
we are saying is that if a particle travels from some localized source to a localized detector by scattering 
twice, first off impurity A and then off B, its path length will be different than if it scatters first off B 
and then off A. As long as this difference in path length is large compared to the spatial extent of the 
wave packet, no significant interference can occur. The basic criterion which is relevant here, and in all 
applications of kinetic theory, is that the mean time between collisions must be large compared to both the 
temporal duration of the initial wave packet (or Ti over its energy), and to the time duration of a single 
scattering event (as determined by the energy derivative of the scattering amplitude). 
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III. LADDER DIAGRAMS AS INTEGRAL EQUATIONS 



A. Relativistic Schrodinger equation 

A concise and convenient method for representing the propagation of nearly on-shell 
particles in a soft gauge background is with the relativistic one-particle Schrodinger equation. 
The free equation is 

id t iP = E^, (3.1) 

where E p = y/p 2 + m 2 and m is the appropriate mass (in our case m^, which characterizes 
the in-medium dispersion correction of a fermion of negligible explicit mass). The corre- 
sponding retarded free propagator is just 

and satisfies 

6{t-t') (pI^MOIp') = f^r e~ ip0 ^ G (P) <p|p'> . (3.3) 

J ATI 

We will use the traditional normalization of one-particle quantum mechanics, 

<p|p'> = (27r) 3 5( 3 )(p-p'), (3.4) 

as opposed to the usual convention in relativistic field theory, which has an extra factor 
of 2E p on the right hand side. Except for the change in normalization, Gq(P) is just the 
positive-energy pole piece of the usual propagators of relativistic field theory (or, equiva- 
lently, the energy denominator of time-ordered perturbation theory when only single particle 
intermediate states contribute). The advantage of the description used here is that it does 
not depend on the spin or type of particle; it can be used to economically describe the 
propagation of hard scalars, fermions, or gauge bosons in soft background fields. 

In a soft background gauge field (containing fluctuations on a scale Q ~ gT < T ~ P), 
the single-gluon vertex is just 

igv-A, (3.5) 
up to corrections suppressed by Q/P ~ g, where v is defined by 

« = £ = (l,v)~(l,$). (3.6) 

A simple, quick way to see this is to replace P by P — gA in the Schrodinger equation ( |3.1| ) 
and propagator (|3.2|) : 

id t ip - gA = E p _ gA ip. (3.7) 

Because the background field is soft, one may treat A and P as if they commute, since 
commutators will replace the operator P, whose relevant matrix elements are 0(T), by 
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Q ~ gT, and hence only generate a subleading contribution. Expanding the equation to 
first order in A, one then obtains 

id t ip = (E p -gvA)ij + 0(A 2 ), (3.8) 

which gives the interaction G3.5|) . Only the "convective" contribution of the charged particle 
to the gauge current appears in the soft gluon vertex ( |3.5| ); the spin dependent contributions 
which would normally distinguish scalars, fermions, or gauge bosons are absent. Spin depen- 
dent contributions to the current are suppressed by an additional power of Q/P ~ g, and 
hence only generate sub-leading corrections to the propagation of a hard excitation through 
the soft background gauge field. 

From the point of view of one-particle quantum mechanics, anti-particles are just another 
type of particle. The propagation of nearly on-shell anti-particles in a soft background field is 
also described by the propagator ( |3.2| ) and vertex ( |3.5| ) except that the gauge field A should 
be in the conjugate representation. Readers who prefer a discussion in terms of standard 
relativistic Feynman rules may refer to Ref. []25|] . 



B. The integral equation 

The contribution to the leading-order photon emission rate from the sum of the ladder 



diagrams depicted in Fig. [13| may be expressed in terms of the solution to a suitable integral 
equation. Analogous summations of ladder diagrams are, of course, well known in many 
other contexts (such as the Bethe-Salpeter equation for bound states), although the form 
of the resulting integral equation intimately depends on the precise form of the propagators 
making up the rungs and side rails of the ladder diagrams, as well as the specific vertex 
factors at the ends of the diagrams. 



The group structure of the diagrams of Fig. [13] is trivial: for every explicit power of 
g 2 , there is also one factor of the quadratic Casimir Cr for the color representation R of 
the emitter. This Casimir is defined by T^T^ = Cr where Tjj are the generators of the 
representation, normalized so that trl^T^ = \ 5 ab . For conventional quarks or anti-quarks 
in the fundamental representation of color SU(3), this Casimir is 

C R = \ [fundamental SU(3)] . (3.9) 

Our first step will be to resum the self-energy insertions of Fig. [TJ] by switching from the 
propagator (|3.2| ) to the resummed propagator 

(3.10) 



P °-E p ±ir l 



where the + sign gives the retarded propagator and the — sign the advanced one. The width 
T p is proportional to the imaginary part of the self-energy of the hard particle. The explicit 
form for this width is discussed below. The corresponding real part of the self energy is 
absorbed into E p . In fact the real part, unlike the imaginary part, is dominated by hard 
rather than soft collisions and cannot be directly computed in the soft background formalism 
we are using. Instead, it can be taken directly from the well-known hard thermal loop result, 
and just gives the usual effective thermal mass. For hard quarks, the effective energy is then 
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FIG. 15. The leading-order contributions to the current-current correlator expressed in terms of 
the solution to a linear integral equation. 



E r 



p 2 + m 2 QO , 



(3.11) 



up to yet higher-order corrections, where [|35|.[36[ 

m 2 0O = \C R g 2 T 2 . (3.12) 
The sum of interferences represented by the diagrams of Fig. [H| may now be rewritten in 



terms of a linear integral equation as depicted in Fig. [15]. To write this explicitly, it is conve- 
nient to work in frequency and momentum space. Consider the contribution to the correlator 
W^D from the x° < portion of the x° integration in Eq. ( [2. Q\ )!^\ If we rewrite the time 



evolution matrix elements in Eq. ( p.7|) as (((p^|Z7(0, x )|p/) (pi\U(0, x ) |Pi)*)) then we may 
use retarded propagators in evaluating these matrix elements. The resulting contribution to 
the current correlator, illustrated in Fig. [15], is 



n(p+k) [1 ± n(p)] (p|^(k)|p+k>* S U {P- K), 



(3.13) 



where S U (P; K), which represents the shaded vertex in Fig. [TR satisfies the integral equa- 
tionQ 



p+k 



P+k, 



(p° + A; ) 

(p\Uk)\p+k) + g 2 C R J (\v p+k -A(Q)] [v p -A(Q)}*js u (P-Q;K) 



(3.14) 



The color matrices and the momentum-conserving; delta function (2vr) 4 5 (4) (- ■ ■) are to be 
understood as already factored out of the ((AA)) correlation in this (and subsequent) for- 
mulae. The x° > region of the x° integration in (|2.6|) simply gives the complex conjugate 



10 We could equally well focus on the a; > portion. As noted below, the two contributions are just 
complex conjugates of each other. It happens to be the x° < portion which leads to the specific form of 
the integral equation derived in our earlier work J25j. 

11 6>„ here is analogous to the combination TT> V of Ref. |25[ . We have chosen our conventions for defining 
S v so that later equations will closely resemble those of Rets. [ p5|2q 1 . 
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(the replacement of retarded by advanced propagators), so that the total bremsstrahlung 
contribution, from a single carrier type and spin state, is 



W *l em (K) = 2 Re 



n(p+k) [1 ± n(p)] (p|j M (k)ip+k)* S„(P; K) 



(3.15) 



Note that the photon four-momentum K is fixed, but the integral equation relates 
S U (P; K) at different values of the four-vector P. In fact, at the order of interest, the 
dependence on P may be simplified. Let p\\ represent the component of p in the direction k 
of the photon, and let p^ be the part of p perpendicular to k. In our process, the momentum 
transfer is O(gT), and so the relevant values of p^ are O(gT). Expanding in both p^ and 
the thermal mass moo, one has 



2 2 

p° - E p + zr p ~ (p° - ^ + m - + *r p , (3.16) 

2 P\\ 

(pO + A: ) - Ep+k + *r p+k ~ (p° -p„) - + ^ P+k , (3.17) 

where we have assumed that the photon is on-shell, and that p\\ > so that the photon is 
traveling in the same direction as the emitter. The last two terms in both equations are 
0(g 2 T), and so each of the propagators will be suppressed unless p°—p\\ is also 0(g 2 T).Q If 
S U (P; K) is only ever going to be integrated against functions that are smooth on the scale 
of 0(g 2 T) variations in p°—p\\, then one may approximate 

( i \* i 2ni 5(p°— py) 9{jp\\) 



(p° + k°) - E p+k + zr p+k ; p°-e p + it p (E p+k + zr p+k ) - (p + e p - ir p 

2n5(p°-p ll )9(p ll ) 



i5E + (T p+]l + r p ) 

where 



(3.18) 



5E es ^4nrr " k ° + E » - E ^ ■ ^ 

2p\\ {k+p\\) 

This is commonly known as the pinching pole approximation, and in our analysis is justified 
for computing the rate at leading order in g. Making this substitution into the integral 
equation ( |3.14j ), one sees that the resulting solution will have the form 

S U (P; K) = S u (p; K) 2tt 5(p°- Pl \) 0(p|,) . (3.20) 

Inserting this form, and rewriting 

S(P°-P\\) 5((p°-q°) - (pH-gy)) = S(P°-P\\) 6(q°-q\\) , (3.21) 



12 If pu < 0, then the two propagators are peaked near p° — —pu and p° + k = \pu + k\, and it is impossible 
for both propagators to simultaneously be (nearly) on-shell. Hence, only the region pii > contributes to 
the leading-order bremsstrahlung rate. 
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one finds 

W^{K) = 2Re / n(p+k) [1 ± n(p)] (p|j M (k)|p+k)* S„(p;k) %„), (3.22) 



where 

S„( P ;k) 1 



(p|j„(k)|p+k) 

V (3.23) 



ijs + (r p+k + r p 

+ < 7 2 C ii | 27r%°-gn)(K +k -A(Q)] [v p • A(Q)]*^ S^p-q; k) 



In terms of the current matrix elements appearing in the original definition fl2.6|) , the function 
S^pjk) is given by that portion of the time-evolved off-diagonal matrix element of the 
current, evaluated in the fluctuating background gauge field, which is phase- coherent with 
the emitted photon. That is, to leading order in g s , 

J d 4 x6(-x°) e~ iKx ((p|^(x)|p+k))) = S„(p;k)0(p||) . (3.24) 

Time evolution is, of course, implicit in j u (x) = j u (x°,x) = U(0,x°) j u (x) U(x°,Q). 
In the interaction term of Eq. ( |3.23| ) one may, at leading order, replace v p • A by 

A + = A Q -A\\ (3.25) 

and rewrite the integral equation in the form 

(p|j„(k)| P +k> = [z^ + (r p+k + r p )]^( P; k) 

-g 2 C R I 2n5(q°- m ) (A + (Q)[A + (Q)r) S„(p-q;k) . (3.26) 

As it stands, this integral equation mixes together different values of p\\ (as well as different 
p_i_) in the argument of 5(p;k), due to the O(gT) momentum transfers in the interaction 
term. The characteristic size of p\\ for a hard quasiparticle is 0(T), and no elements in the 
equation are sensitive (at leading order) to gT variations in p» . Consequently, we can treat 
p\\ as fixed inside the interaction term,[^] replacing the equation by 

<p|j„(k)| P +k> = [z^ + (r p+k + r p )]^( P± ;p||,k) 

-g 2 C R I 27T6(q -q ll )(A + (Q)[A + (Q)r}S I/ (p ± -c l± ;p h k). (3.27) 



13 This may be argued in more deta il as follows. Given the linearity of the integral equation, one could 
multiply the source ( p|jV( k)|p+k) in ( 3.26 ) by S(p\\— P\\'), at the cost of introducing an integral over p\\' in 
the final expression ( |3 ■ 22| ) . The resulting solution to the integral equation would then have its dominant 
support on values of p\\ which differ from pu' only by total momentum transfers of O(gT). Because the 
(p|j„(k)|p+k)* 9{p\\ ) factor in (3.22) is smooth in pn (for pii of order T), no error is made at leading order if 
P|| is replaced by p^' in these terms. If the modified ( |3.26 ) is then integrated over py, and J dp\\ S(p;p\\', k) 
is renamed S(p±;pu' ,k), one arrives at Eqs. (3.22) and (3.27). 
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This is now a linear integral equation in p± space, for fixed values of p» and k. 

One remaining complication is that both the width T p and the total scattering cross sec- 
tion have infrared divergences in our approximations, proportional to the Q — > divergence 
of Jq 27i5(q°— q\\) ((A + (Q)[A + (Q)]*)) ■ This divergence of T p in HTL perturbation theory is 
well known ]37| , arises from collisions via the exchange of unscreened low- frequency magnetic 
gluons, and is a symptom of sensitivity to ultra-soft non-perturbative magnetic physics at 
the momentum scale g 2 T. But this apparent infrared sensitivity is actually illusory in our 
problem, as may be seen by rewriting the width in a form similar to the interaction term in 
Eq. ( 3.27[ ). The width comes from the imaginary parts of the self-energy insertions in Fig. 



p o=E p 



|M| for nearly on-shell particles and, to leading order, is given by 

T P = to ]pc R J Q (AHQ) A-(-Q)] 

- </C R j rcS( q °- q ,) (A+(Q) [A+(Q)]'). (3.28) 
Substituting this into ( |3.27| ), one obtains the infrared-safe equation 
(p\j u (k)\p+k) =z5ES„(p ± ; P \\,k)+g 2 C R [ 27r5(q - qil ) (A + (Q) [A + (Q)Y 



x [Sv(p±;p\\,1z) - S , „(pj_-qi;p||,k)] . (3.29) 

We may also tidy up the expression for W^ u by noting that n(p) ~ n{p\\) [since p_i_ is O(gT)], 
so that to leading order 

W^ m (K) = 2Re / n(p l{ +k) [l±n(p {l )\ (p|j M (k) |p+k>* 5,(p ±;P ||,k) 0(p\\) . (3.30) 



C. Current matrix elements 

The current matrix element (p|j M (k)|p+k) is the only ingredient of our equations that 
depends on the nature of the emitter (i.e., scalar or fermion). Since p^ is parametrically 
small compared to p\\ , it is sufficient to evaluate these matrix elements in the limit of small 
transverse momentum. We only need the transverse components of because the currents 



are dotted into photon polarization vectors in Eq. ( |2.2|) . By rotational invariance about the 
k axis, these must be proportional to p_i_. Let qe denote the electric charge of the emitting 
particle (so that up quarks, for example, have q = 2/3). Then in the small p^ limit, the 
transverse current matrix element takes the form 

(p|j±(k)|p+k) ~ 2qep ± J p (3.31) 

where we have chosen to factor out the charge of the emitter. The explicit form of the 
"splitting function" J p ^ p ^ + k, which depends on p\\ and k as well as the spin of the emitter, 
will be discussed momentarily. Given the linearity of the integral equation (|3.29| ) (and the 
fact that it only couples different values of p_i_, not p\\), we can factor out all dependence on 
the emitter type by redefining the transverse part of S u (p±;p\\, k) as 
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Sj_(p±;P||,k) = qej p ^ n+k f(p ± ;pn,k) 

to obtain 

E«<rW = 2(ge) 2 Re 



n( 



-k) [l±n( P]l )]\J. p +& | 2 2p ± -f( Pj 



(3.32) 

■ k)*(Pll), 
(3.33) 



where f (pj_;p||, k) is the solution toQ 

2p ± = z^f(p ± ; P ||,k)+/C i? / 27r5(g°-g||)((A + (g) [A+(Q)]* 

X [f(p±;p||,k) -f(p ± -q ± ;p||,k) 



(3.34) 



This gives the result for bremsstrahlung production from a single type and spin (and initial 
color) of charged particle in a way that clearly isolates (i) the dependence on the type of 
the particle in the splitting factor | t ^,..«_ p .. + fc| 2 , and (ii) the dependence on the details of the 
frequency-dependent correlation of the background field in the correlator ((AA)). 
If the charge carriers are scalars, one would have 

(2p + k)_L 



(p|j±(k)| P +k) = qe 



(2E p )y^2E p+k )y^- qe \ P]l (p l] +k)}^ 

The numerator in the middle expression is the usual photon vertex in relativistic normal- 
ization, while the denominator arises from our use of a non-relativistic normalization ( |3.4|) 
for our states. Making a similar calculation for fermions, one obtains^ 

1 



(3.35) 



\Jt 



Pii<-Pu+fc| 



4p|| (p||+fc)' 
P\\ 2 + (P\\+k) 2 



scalars; 



fermions. 



(3.36) 



8pii 2 (pii+£0 2 

This result can also be expressed in terms of the leading-order Altarelli-Parisi (or DGLAP) 
splitting functions for (hard) photon bremsstrahlung from a charged particle q, 



\Jt 



Pll+-PH+fe| 



kP q(l) ^ q {z) 



l (p\\+k) 



(3.37) 



=Pn/(Pn+fe) 



whereQ 



14 Rotation invariance about k implies that f(p_i_;pii, k) m ust eq ual p^ times a scalar function of |p_i_|, p\\ 



and k. (See Ref. |26|| .) But the resulting integral equation ( p.34|) is more compact if f(px;pii,k) is left as a 
vector function. 

15 See, for example, Ref. pa] for more detail, bearing in mind the change to the non-relativistic normal- 
ization used here. 

16 -^f(s)- 1 — 1 = + z2 )/(l — z )+ + | <5(1 — z )\ i s the standard result for gluon bremsstrahlung from a 
quark. The corresponding splitting function for photon production is obtained by replacing Cr by 1. The 
treatment of the soft photon singularities at z=l is not relevant to us, since our photons are hard; hence we 
may ignore the 5(1— z) and the endpoint prescription on the (1— z) denominator. 
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2Z scalars; 



1-z' 



fermions, 



z 



for z < 1. 



D. Relationship to Migdal's equation 

Let us note in passing the relationship of our result to similar equations written by 



Migdal P9| , p0| for the case of hard bremsstrahlung during electron scattering from the static 



Coulomb fields of randomly distributed atoms. If the background were static, we could write 
[A + (Q) [A + (Q)T} = 2tt%°) (<V(q) L4+(q)]*}, (3.39) 
and the collision term in (|3.29|) would then become 



g 2 C R / q 27r%||) (<V(q) L4 + (q)]*)) [S u (p ± ; P h V - ^(p ± -q ± ; p h k)] . (3.40) 

Migdal has an equivalent expression (Eq. (9) of Ref. p0| ) but with 8(q\\) replaced by 

1 6(E p+tl - E p ) + \ 5(E p+k+Ci - £ p+k ). (3.41) 

At leading order, however, both of these delta functions are equivalent to 5(q\\), and so our 
result reduces to Migdal's in this special case. 

Migdal went on to solve his equation in the approximation that the LPM effect was 
parametrically large — roughly, that SE <C T p (or equivalently that the formation time was 
large compared to the mean free collision time), in which case it turned out he was able to 
expand in the logarithm of the ratiof£] \n(T p /8E). In our application, 5E and T p are of the 
same parametric order, and such an expansion is not useful. 

IV. PAIR ANNIHILATION 

The treatment of pair annihilation is very closely related to that of bremsstrahlung — a 
fact which is slightly obscured in the present time-ordered discussion but is more manifest in 
our original treatment of Ref. . The contribution to the current-current correlator 
due to pair annihilation, in a form analogous to Eq. (EOjp for bremsstrahlung, is 



Wp h (K) = ij d 4 xe-' & J^n(p)n(p) (pp|j M (0)|vac) (vac|j^)|pp>\\ . (4.1) 



17 At very high energies, this logarithm is modified in Midgal's analysis by the finite nuclear size of the 
atomic scatterer, which is also irrelevant to our application. 
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t=0 



t-x 



t=0 



t-x 




FIG. 16. Two (time-ordered) diagrams whose interference contributes to the pair annihilation rate. 

Using the fact that particles propagate independently in our approximation (the only inter- 
actions are with the background field), this can be written in terms of one-particle evolution 
amplitudes as 



d x e 



-iKx 



n{p) n{p) (pp|j M (0)|vac) (vac|^(x)|p'p') (p'p'|[/(x u , 0)|pp) 



d 4 x e 



iKx 



pp'pp' 



nip) nip) (pp|j M (0)|vac) (vac|j„(x)|p'p') 

x ((p'|C/(x°,0)|p> (p'|C/(x°,0)|p)\. 



Translation invariance now gives the momentum conservation constraint 

k = P + P = p' + p'. 



(4.2) 



(4.3) 



Note that we have ignored processes where the particle and antiparticle interact directly 
with each other via gluon exchange, rather than with the random background field created 
by other color charge carriers in the plasma. Such processes are sub-leading in g, as are 
vertex corrections to bremsstrahlung.[^] 

The pair annihilation contribution to the current correlator W^ v represents the interfer- 



ence of the two processes shown in Fig. [16]. The interference of these amplitudes can again 
be depicted by a diagram like Fig. |^, but now the interpretation is slightly different. The left 
photon vertex represents the conjugate of the first diagram of Fig. [H], and everything else 



represents the second diagram of Fig. 16 . The same time and space ordering considerations 



18 A simple qualitative way to understand this is to consider the typical interaction energy times the 
photon formation time (which, among other things, is the time the annihilating pair spend within a Debye 
screening length of each other) in the center of mass frame of the almost-collincar annihilating pair. If this 
product is small compared to 1, then the effect of direct gluon interactions between the pair is suppressed. 
The formation time is \jg 2 T in the plasma frame, and becomes smaller by a factor of m^jT ~ g when 
boosted to the center-of-mass frame, giving 1/gT. The characteristic distance between the pair over the 
formation time is r ~ 1/gT, giving a characteristic potential energy of order g 2 /r ~ g 3 T. The product is 
then order g 2 -C 1. For comparison, the interaction time and distance for interaction with a random colored 
particle are both order 1/gT, and the energy-time product is again g 2 . But the pair encounters N ~ l/<? 4 
such random particles during the formation time [number density T 3 times Debye screening volume 1/(<?T) 3 
times formation time \jg 2 T divided by time 1/gT interacting with each particle]. These encounters will 
on average cancel each other in their effects except for statistical fluctuations, which will enhance our g 2 
contribution per random encounter by y/~N to give a net, unsuppressed effect of g 2 ^/N ~ 1. 
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discussed for bremsstrahlung again imply that the only relevant correlations at leading order 



are those depicted by Fig. |T3[ The resulting integral equation is closely related to the one 
for bremsstrahlung: 



W?f{K) = 2Re / n(k-p)7i(p) (vac|j M (k)|p, k-p)* S„(P; K) 



(4.4) 



withQ 

S U (P;K) 



(k -p°)-E^ p +iT k _ p J \ p 0-E p + iT p J 
(p,k-p|^(k)|vac) + <7 2 C* [ f[v^ p -A(Q)\ [v p ■ A(-Q)]\ S V {P-Q;K) 



(4.5) 



The pinching pole approximation is now 



(fc°-p°)-£k- P + TV f 



p°-e p + ir. 



where 



5E 



P± + ™>lo 



(k- 



2iri5(p°-p\\) 9{p\\)9(k-p\ 
(£ k _ p + ziV p ) + {E p + tT p ) 

27rS(p°- Pll ) 6( Pll ) 6(k- Pll ) 

iSE + (r k _ p + r p ) 



-A; - £ k _ p - E p . 



- k° 



(4.6) 



(4.7) 



The same simplifications as for bremsstrahlung then reduce these expressions to 
E« W^{K) = 2 (ge) 2 Re / n(k- Pl] ) [1 ± n(p„)] \J„ hk - n ? 

x 2p± -f(p±;p||,k) e(p\\)6(k-p\\), 

and 

2 P± = t5E{(p ± ;Pii,k)+g 2 C R [ 2vr 5(q°- qil )(A + (Q) [A+(Q)} 

J Q 

x f(p±;p||,k) -f(pj_-qj_;pn,k) 

where we have defined the small p_i_ limit 

(vac|j M (k)| P) P — k) ~ 2ge pj_ jTvac^-pii^-pn 
analogous to Eq. ( |3.31| ). These "joining functions" are given by 



(4.8) 



(4.9) 



(4.10) 



19 A note on signs: the i's from the two soft gluon interactions cancel the minus sign associated with the 
fact that the color generators for the anti-particle are — T* instead of T. 
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4p|| (k-p\\)' 
" P| " fc " P " r " | P|| 2 + (A:-P||) 2 
- 8[p|| (^-P||)] 2 ' 



(4.11) 



scalars; 
fermions; 

— p\\ in the bremsstrahlung formula 



which is the same, up to an overall sign, as taking pn — > - 

O). 

The results (|4.S| ) and ( |4.9| ) can be cast into exactly the same form as the bremsstrahlung 
results ( |3.33|) and ( |3.34|) , with the only difference being a factor of 9(—p\\) 6{k + p\\) in Eq. 
( fi~8l) instead of 0(p\\). To see this, everywhere redefine p — > — p in Eq. (|4.8|) , including 
P\\ ~ * ~P\\- Then define 



f fa±,Pf k ) = f (-P.L, ~Pf k )- 

The equivalence is completed by the identity 

n (~P\\) = =F[1 ±™( 



(4.12) 



(4.13) 



and the relationship between the squares of current matrix elements ( 3.36Q and ( 4.11 ) just 
noted. 



V. FINAL FORMULAS FOR PHOTON EMISSION 



To produce a final formula for the leading order contribution to hard photon production 
from near collinear bremsstrahlung and pair annihilation, we just need to combine the 

for pair 



defining relation (2.1) with the result (|3.34j ) for bremsstrahlung and the result 
annihilation, transformed as just discussed, and finally sum over the colors, spins, and flavors 
of possible emitters. The final result can be neatly packaged in the form 



dv L ™ 
d 3 k 



d 2 P± 

(2vr) 2 



Uf(p\\+k) [l-7lf(p||)] 

Ref(pj.;pii,k), 



-fc| 2 2 PJ 



(5.1) 



where f(pj_;p||,k) satisfies the integral equation (|3.34p , repeated here for convenience 



2 P± = i5Ei( P±] p h k)+g 2 C R f 2n5(q - qil )(A+(Q) [A + (Q)]* 

J Q 

x f(p±;pi|,k) - f(p ± -q ± ;p||,k) 



(5.2) 



and SE is defined in Eq. ( ^.19| ). In the result ( |5.1| ), the sum over s is a sum over the 
relevant species (flavors) of quarks (i.e., u, d, s, c, . . . ), and the current matrix elements are 
given by the fermionic case of Eq. (|3.36|) . N s = 4 counts the spin and particle/anti-particle 
states for each flavor of quark, and da ; = 3 is the number of colors of a quark. As usual, 
ri{(ou) = \j[eP^ + 1) is the Fermi distribution function. In the result fl5.1|) , the contribution 
to the p|| integral from p\\ > corresponds to bremsstrahlung from quarks, the p\\ < —k 
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contribution is identical (after substituting py — > —A; — P||) and accounts for bremsstrahlung 
from anti-quarks, and the —A; < p« < contribution corresponds to pair annihilation. 

All dependence on the nature of the scatterers which influence the photon emission 
rate is contained in the relevant ((AA)) correlator. As discussed earlier [c.f. Eq. ( [2.4|) ], the 
background field correlator is given by 2[nt,(q°) + 1] ImCrR, e t(Q), with the gluon propagator 
evaluated in the hard thermal loop approximation for soft momenta QcT. Using standard 
results for the hard thermal loop self-energy, one finds the gauge-invariantFl correlation 



A + (Q) [A + (Q)Y 




(u/qY 



(9°) 2 + n T (Q)F 



(5.3) 



when q° is set equal to q\\ by the delta-function in Eq. ( [5.21 ). Here, mrj is the (lowest-order) 
color Debye mass, given by m 2 D = (1 + |iVf) g 2 T 2 for QCD with Nf flavors, ancQ 
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Together with Eqs. (|5.1|) and (|5.2[), this reproduces our previous results of Ref. p5 



VI. GLUON EMISSION 

We now wish to extend of the above analysis to the case of gluon emission. Other 
authors who have considered the LPM effect have also discussed gluon emission f20|-|24 



Our generalization is similar to these previous treatments, but our final result is the first to 
give the full correct leading-order answer at high temperature because it is the first to treat 
correctly the dynamics of the scatterers. 

The key difference from photons is that the gluon also carries color. So, whereas the 
photon merely contributed a phase exp(ik°t) to the time evolution amplitude, an emitted 
hard gluon will feel the random colored background field, just as the emitting particle does. 
But provided we consider hard gluon emission, the emitted gluon is distinct from the soft 
background, and should be treated as another (quasi)particle. The analog of Fig. |7| for 
gluon emission is shown in Fig. [17]. The analogs of Fig. [5] and Fig. [13] are shown in Fig. [L8|. 
Much of the argument for photon emission still holds; the gluon and quark must be nearly 



20 The change in the (A + (Q)A + {Q)*) correlator under an Abelian gauge transformation is proportional to 
q + = q°— qu, but we are evaluating at q + — 0. Hence this change vanishes for any transformation whose 
Fourier transform is non-singular at q + = (e.g., any transformation of compact support). Because gauge 
field fluctuations on the gT scale are perturbative, the non- Abelian part of a gauge transformation changes 
the correlator by a subleading amount suppressed by g. 

21 In the convention used here, the hard thermal loop self-energy n A " y ((5) = (Q 2 /q 2 ) ITl(Q) + 
^t{Q)'P^ u , where "Pl and Vt are longitudinal and transverse projection operators. Various other authors 
(sometimes including ourselves) instead use the notation IIl(Q) to denote our — [Q 2 /q 2 ) IIl(Q). 
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t=0 



t=x 




FIG. 17. Two diagrams, time ordered from left to right, whose interference contributes to the rate of 
gluon bremsstrahlung. The only difference from the photon emission case of Fig. is that the gluon, as well 
as the emitter, can interact with the soft background field. 





FIG. 18. A single diagram depicting the interference of the two gluon emission amplitudes of Fig. 17. 
The left-hand version depicts the product of amplitudes in a given background field; the right-hand version 
shows the result after averaging over random background fields. Individual background field correlators can 
connect any two of the three hard lines, and double lines represent the rcsummcd propagators of Fig. nil 



collinear, and consequently the interactions are ordered in the same way as before. The 
new complications are, first, that there is a color matrix at the hard particle vertex, and 
second, that there are now correlations between the soft gauge field felt by the gluon and 
either emitter line, not just between the two emitter lines. 

Because the interactions are ordered as before, it is still possible to resum the diagrams 
by an integral equation similar to Eq. ( |3.l4| ). The only difference is that there are now three 
elements in the collision term, corresponding to the three kinds of correlations between lines 
shown in Fig. 



18. The group theoretic coefficients are easily found with the help of 
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where T R denote representation R color generators. Since the soft gluon correlators are 
ordered in time, this is sufficient to determine the group factor for the whole diagram; each 
line between emitters gives a factor of (Cr — |C*a) and each line from an emitter to the 
emitted gluon gives \Cp^. 

The other complication is that, whereas before a cross- rung always changed p and left 
k the same, now it can either change p, change k, or change both. This is a reflection of 
the fact that the gluon, unlike the photon, can scatter during the 1 / g 2 T time scale of the 
process. In the photon case, we had a natural fixed direction k with respect to which we 
defined (whose changes we had to keep track of) and m (whose changes we could ignore). 

One inelegant but concrete possibility would be to proceed as before but (i) pick by some 
convention any direction nearly collinear with k, p and p + k to define the "parallel" (||) 
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direction, and (ii) keep track of changes to k^ as well as p_i_. For example, consider picking 
the || direction to be that of the momentum of the emitted gluon line at the right-most 
corner of Fig. 18. The analog to our results for photon emission would then be0 
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where we have included a radiation stimulation factor l+n b (k), defined A pk = p_|_ — p||kj_/fe|| 
as the component of p perpendicular to k (at leading order), and where the accompanying 
integral equation is 
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(6.3) 



The three terms in the integral correspond respectively to gluon exchange between the top 
and bottom rails, the top and middle rails, and the bottom and middle rails of Fig. 18). The 
momenta flowing from left to right in the top, middle, and bottom rails of that figure have 
been labeled p, k, and — p — k, respectively. The energy difference SE (for which we have 
suppressed the arguments pj_, kj_, p|| and k\\) is 
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(6.4) 



This expression incorporates the fact that in an equilibrated thermal plasma, the hard gluon 
dispersion relation includes an effective mass m 2 = m^/2. 

The preceding is a valid but inelegant approach because it does not exploit the rotational 
invariance of the problem. Other than |f s |, the elements of the integrand in the rate formula 
( |6.2| ) are invariant under small angle (6 ~ g) rotations, which do not change m and kn at 
leading order, but do change p^ and kj_. So one could split the d 2 p± d 2 k± into one integral 
with respect to the relative momentum A pq , which appears explicitly in ( |6.2| ), and another 



22 <fTg PM /d 3 fc represents the differential rate per unit volume of processes that increase by one the number 
of gluons with momentum k. One might be tempted to think that the integral of this quantity is therefore 
the rate per unit volume of processes that increase the total gluon number by one. This is incorrect since 
pair annihilation of hard gluons contributes to dTg Pm /d 3 k but decreases rather than increases the total 
number of gluons. In addition, once one integrates over k, one has to avoid double counting final states 
in the bremsstrahlung process. We are unaware of any simple physical interpretation of the integral of our 
(ffg PM /d 3 k, although this differential rate is closely related to ingredients needed in the evaluation of various 
physical quantities, such as leading-order transport coefficients jlTj]. 
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integral over small rotations. The integration over small rotations could then be absorbed 
into f s , defining a new function that depends only on A pq . It turns out that the most 
symmetric procedure does not use A pq but rather the quantity h = fei|A pq , which can be 
written as 

h = A;|| pj_ — _p|| k_L. (6.5) 

This may also be written as h = (k x p) x en, where eii is the unit vector in the || direction. 
The corresponding equations that only track changes in h are 

dT LPM a s AT , f+°° dp f d 2 h . , , ri . . . ri 

-sfc- = is? ? N - d - c - L i I w? n - (p+k) 11 T nM] 11 + " b{k)] 

x Uiir+kf 2h ' ReF,(h;p, k) , (6.6a) 

where F s (h;p, k) is the solution to the integral equation 

2h = % 5E(h; p, k) F s (h; p, k) + g 2 [ 2n 5(q°- q]l ) (A + (Q) [A+(Q)}*} 

x |(C S - |C A ) [F s (h;p,k)-F s (h-kq ± ;p,k)} 

+ \C K [F s (h; p, k) - F s (h+(k+p)q ±] p, k)} 

+ \C A [F s (h; p, k) - F s {h-p <u;p, fc)]|, (6.6b) 

with 

777^ 777^ 2 

<B (^*, = 1 . + ^ + 5F -^- SEF _ j5 . (6.6c) 

To simplify the appearance of these results, and the equations below, we have dropped the 
|| subscripts on p\\ and k\\. It should be emphasized that evaluating these expressions does 
not, in fact, require selecting any particular convention for defining the || direction. Also, 
note that rotation invariance implies that F s (h;p, k) must equal h times a scalar function 
of |h|, p\\, and k, but the integral equation (|6.6b| ) is more compact if F s (h;p, k) is left as a 
vector function. 

In the result ( |3.6a| ), the sum over species s now includes gluons, which have N s = 2, 
d s = g?a, and C s = Ca (with g?a=8 and Ca=3 for QCD). The splitting function for gluons is 

1 + ** + (!-*)' (67) 



for z < 1, so that 



l7 ( 9 ) l2 _ P 4 + k* + (p+ky 



This expression, when multiplied by k 3 , is symmetric under interchanges of the three gluon 
momenta p, k, and —p—k. 
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Further symmetry of our result (whether for quark or gluon emitters) can be made 
manifest by rewriting group factors so that the expression in curly braces in the integral 
equation ( |6.6b| ) is 

\{C R , + C Rl - C R2 )[F s {h;p, k) - F s (h-p 2 q ± ;p, k)] + 
1 1 (C Rl +C R2 -C R3 )[F s (h;p,k)-F s (h-p 3 qr,P,k)} + 
l(C R2 +C R3 - C Rl )[F s (h;p,k) -F s (h- Pl qr,P,k)]y (6.9) 

where (pi,P2,P3) denotes the three momenta (p, k, —p—k) and (C Rl ,C R2 ,C R3 ) denote the 
quadratic Casimirs (C s , C*a, C s ) of the corresponding color representations.^ 

The previous photon emission result can be extracted from the final gluon emission 
expressions ( |6.6| ) by setting Ca=0 inside the integral equation (which corresponds to making 
the emitted particle colorless rather than in the adjoint representation), replacing C s in 
Eq. ( |6.6a| ) with gj?, and using the photon, rather than gluon, thermal mass in the expression 



( |6.6c|) for 5E. (The resulting photon thermal mass contribution to SE is irrelevant for hard 
photons P5fl .) 

Finally, we make a last note about the cancellation of infrared divergences in our result 
for gluon production. The fact that the collision term in ( 3.6b|) only involves differences of 



the form F s (h;p, k) — F s (h — PiQ±; p, k), which cancel in the far-infrared (q_i_ — > 0) limit, 
results from including all three ways of connecting the three rails in Fig. [L8] together with 
resumming self-energies (as in Fig. |TJ]) on all three rails. In contrast, if we had truncated the 



hard gluon line and evaluated the earlier diagrams of Fig. [13] with external color currents, 
these cancellations would have been absent and our result would have had far-infrared diver- 
gences signalling sensitivity to ultra-soft momentum scales. The sensitivity of such truncated 
diagrams to ultra-soft physics is relevant to the evaluation of quantities such as color con- 
ductivity,^ but not to the problem of the hard gluon emission rate at leading order. 
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23 There is a general argument that gives the structure of the group factors shown in Eq. ( jfjty even if i?i , 
R2, and i?3 were any three color representations. Color conservation at either end of Fig. |lq means that 
these three representations have to be combined into a color singlet. Acting on that color singlet state with 
color generators TJ^ + Tg 2 + TJ| 3 should give zero. So acting with Tj^ + Tg 9 has the same effect as acting 
with — Tyj . Squaring both sides and summing over the basis a of the algebra, one then obtains that TfT$ 
(the factor associated with exchanging a gluon between particles 1 and 2), acting on the color singlet, is 



equivalent to multiplication by \{Cz — C\ — C2). The other group factors in (6.9) follow by permutation 



24 For a discussion of color conductivity in t erm s of ladder diagrams, see Ref . pq | . For an equivalent 
discussion in terms of kinetic theory, see Ref. [p9| and the pathbreaking work of Ret. jhj. See also Refs. 
}4l"|,fi"2"f and the original, early work of Ref. [H . 
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